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1. Introduction 

Understanding the scattering of cosmic rays in 
the interplanetary and interstellar plasmas is a 
problem of central importance in astrophysics. The 
diffusion processes of charged particles in the di- 
rections, parallel and perpendicular, to an ordered 
magnetic field (e.g., the magnetic field of the Sun) 
can be described by the diffusion tensor, whose 
components can be (1) calculated using analyti- 
cal transport theories; (2) extracted from numeri- 
cal test-particle simulations; and (3) obtained from 
heliospheric observations. Understanding of such 
observations is a key subject of space physics [1, 

An important way to test analytical theories [_ 
is numerical Monte-Carlo simulations that op- 
erate under the same restrictions imposed on an- 
alytical calculations — stationary or static turbu- 
lence, prescribed turbulence geometry and power 
spectrum, no back-reaction of the particles on the 
turbulence field. Hence, such simulations are called 
"test particle" simulations although, using the same 
approach, other effects such as magnetic field line 
diffusion can be investigated. 

In this Note, the special case of isotropic magne- 
tostatic turbulence will be investigated, which is an 
important test case for both numerical and analyti- 
cal approaches 0, Hol - [l2l | . It will be shown what the 
basic turbulence properties are and how magnetic 
turbulence is usually generated in numerical simu- 
lations (Sec. [5]). Some problems will be discussed 
that are inherent in the basic formulation of numeri- 
cal turbulence, because not all physical requirement 
such as vanishing magnetic divergence and isotropy 
in position and wavenumber space can be fulfilled 
at the same time. Finally, several simulation results 
will be compared and discussed (Sec. [3]). 
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2. Isotropic Turbulence 

Isotropic turbulent magnetic fields can be 
thought of as a superposition of plane waves with 
random phase angles and random orientations. In 
the limit of an infinite number of plane waves, the 
resulting turbulence is homogeneous and isotropic 
[l3[ . There are a number of both analytical and nu- 
merical constraints and also pitfalls. Consider each 
in turn. 



2.1. Analytical constraints 

In general, homogeneous turbulence [l3| is de- 
scribed using a stochastic approach that is based 
on a two-point, two-time correlation tensor 

(Bt(x, t) B* m {x', t')) = R, m (a:, x', t, t'), (1) 

where B\^ m refers to the turbulent magnetic field 
components with l,m £ {x,y,z}. 

To account for a power spectrum that is (at least 
partially) known in wavenumber space 14, [TBI, a 
Fourier transform is applied, where the assumption 
of homogeneity leads to a delta function S(k — k') 
13, IH 12 1- With the additional assumption of a 
time-independent turbulence field (i.e., magneto- 
static turbulence), the result reads 



(B(k) B*(fc')) = S(k - k') P, TO (fc). (2) 

For homogeneous and isotropic turbulence, the 
correlation tensor has the form 13, lj| lj| 
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(3) 



with a(k) G [—1, 1] the magnetic helicity (usually 
assumed to be zero, with some noticeable excep- 
tions) and ei mn the Levi-Civita tensor. 
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The normalization of the correlation tensor P; 



is given through the condition 16 1 



3 

B 2 = Bl + Bl +B* a = d 3 fc P «( fc )' ( 4 ) 

J i=l 

where B corresponds to the (average) turbulent 
magnetic field strength. 

Consider now the three constraints for the 
isotropic turbulent magnetic field. 

Wave vectors. Isotropic turbulence means that 
each orientation of the wave vector has equal prob- 
ability. Consider the Fourier transform (for illus- 
tration purposes in two dimensions) of an isotropic 
function 

F(k) = { d 2 r F(r)( 



= / drrF(r) / d<t> e ikrcos{ -^\ (5) 
Jo Jo 

with polar coordinates k = (kcostp, ksmi/j) and 
r = (r cos 4>, r sin (ft) . Then [l9j 



a ikr cos(i/' — (p) 



n=-oo ^° 

27rJ (fcr) (6) 



(i>-<t>) 



so that the result does not depend on the orienta- 
tion of the wave vector as described through the 
angle ip. Here, J n denotes the Bessel function of 
the first kind of order n. 

For the three-dimensional similar calcula- 

tion is slightly more involved. However, one can al- 
ways transform to a new coordinate system (r, a, 0) 
where a is defined through k ■ r — kr cos a. Then 
/3 is unused in the Fourier integral so that the two- 
dimensional case is recovered. 

Field strength. According to Eq. Q, the mean 
value of each individual field component is deter- 
mined through^ 



Bf = Jd 3 kP u (k). 



(7) 



Due to the fact that the isotropic turbulence ten- 
sor from Eq. ([3]) does not distinguish, for exam- 
ple, the z from the x, y directions it is immediately 



1 Strictly speaking, k = must be excluded because it 
represents a uniform magnetic field. 



clear that, on average, (B 2 ) — {B' 2 ) — (B 2 ). Any 
isotropic turbulence generator must therefore fulfill 
the constraint of equal amplitude field components 
in all three spatial directions. 

Divergence. Every magnetic field must obey 
Maxwell's equation of vanishing divergence, i.e., 
V • B = corresponding to the property that no 
magnetic monopoles exist (although, in the realm of 
quantum effects, the subject remains under active 
investigation [20| V 

2.2. Numerical turbulence generation 

Following the ideas of 5- 8] , the turbulence in the 
Padian code [9( is generated via a summation over 
N plane wave modes as 

JV 

B(x,y,z) ^Re^T^Me^'+H (8) 

n=l 

where (3 n is the phase angle of the plane waves. The 
vector £ n denotes the amplitude direction of each 
wave and is defined as 



= cos{a n )e x >, n +isin(a n )e y /, n . 



(9) 



where a n is the polarization angle. 

The unit vectors e x / jn and e y ^ n are given by 
the first and second lines, respectively, of a three- 
dimensional rotation matrix 



cos tt cos cos V sin < 

Mm = I — Sll) r > ens < 

sin U cos sin sin < 




(10) 



Physically, one can think of the angle a describ- 
ing wave types that vary between fast-mode waves 
(a = 0) and Alfven waves (a = 7r/2), but since no 
time-dependence is considered, the analogy is lim- 
ited. 

The direction of propagation of the plane waves, 
i.e., the z' direction, results from z' = xA^x + 
V ^33- Because the z' direction is always per- 

pendicular to £ n , one immediately has k n ■ £ n = 
for every mode n, which corresponds to V • B = 0, 
thus ensuring that the turbulent magnetic field is 
divergence free. 

For each summand n, all angles 9, 4>, a, and (3 
are randomly generated^ but are then kept fixed. 



2 Note that not the angle 9 but instead its cosine, r\ = 
cos 9, is uniformly distributed. This ensures that the density 
of wave directions is equal for all solid angles d(p drj. 



Thus, the same B results for the same set of coor- 
dinates (x,y,z), corresponding to what is called a 
"turbulence realization" . 

By constraining the angles 6 and a @, [§] , non- 
isotropic turbulence geometries such as slab and 2D 
can be obtained. To fulfill the additional constraint 



that B _L e z [2l|, |22j , the polarization angle has to 
be set to a = it/2. 

The amplitude function A(k n ) is defined through 



N 



A 2 (k n ) = G(k n )Ak n J2 G(ku)Ak u 



(11) 



where, for example, the turbulence spectrum G(k n ) 
is of the form 1231 



G(k n 



(i + K?) 



(s+q)/2 ■ 



(12) 



where q and s are the energy range and inertial 
range spectral indices, respectively. A logarithmic 
spacing of the wavenumbers is commonly used so 
that Ak n /k n is constant. 

2.3. Normalization of the turbulence 

The turbulent field as generated through Eq. © 
should be normalized to unitylf] Therefore, it is 
required that, on average, 



B 2 = (Bl + Bl+Bl) = l. 



(13) 



If such were not the case, the resulting transport 
parameters would be falsified, i.e., would be too 
large (small) if the magnetic field strength were to 
be smaller (greater) than unity. 

However, on implementing Eq. ([5J as is, one finds 
that: (i) the average strength of the turbulent fields 
is considerably smaller than unity; (ii) the B z com- 
ponent is (also on average) smaller than the other 
two components, which contradicts the requirement 
of isotropy. Both of these drawbacks are due to the 
fact that 



( 



cos a cos tt cos < 



cos a cos tt sin < 



i sin a sin < 



I sin a cos ( 
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2 ) = 1/4- 



3/8 
3/8 
(14) 



3 Technically, such is due to the requirement that a unit 
magnetic field vector eg is generated, which, in the equation 
of motion, is "manually" scaled with the factor B/Bq to 
yield the requested turbulence strength relative to the mean 
magnetic field strength, Bq @. 
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Figure 1: (Color online) The parallel scattering mean free 
path as a function of the particle speed. The errorbars with 
red triangles and blue dots show the simulation results from 
the Padian code [tj with and without renormalized wave 
vectors, respectively. Previous simulation results [(J do not 
use any normalization of the turbulent magnetic field (black 
squares). The black dashed line shows an analytical result 



If the turbulent magnetic field components are di- 
vided by the mean values of the unit vector compo- 
nents, i.e., by \/3/8 and by 1/2, respectively, then 
all components of the turbulent magnetic field have 
equal means and the total magnetic field strength 
is approximately unity, as required. 

2.4- Fulfilling physical constraints 

One is faced with a choice because one has three 
options: 

• using the original turbulence generation mech- 
anism; but the original form for the turbulence 
is not in agreement with equal mean values 
for the three magnetic field components, as re- 
quired for isotropy; 

• normalizing the turbulence as described above; 
but the normalization factors are not compat- 
ible with the requirement that the divergence 
of the turbulent magnetic field be zero (note 
that one still has (V • B) = on average); 

• renormalizing the turbulent wave vector via 
k x ,y — I k Xi yy/3/8 and k z — > k z /2 so that 
k n ■ = Vn is restored; but then the wave 
vector is not isotropic any more. 

The important point to note is that, using the 
approach described by Eq. (O, not all three re- 
quirements can be fulfilled at the same time. Such 



is comparable, e.g., to various formulations of 
smoothed particle hydrodynamics 25|, [26|, where 
conditions such as V • B — or the conservation of 
mass, energy, or angular momentum are frequently 
violated. It is therefore left to the resulting trans- 
port parameters to decide which option gives the 
best results. 



3. Results and Conclusion 

In Fig. [TJ the parallel mean free path is shown 
as resulting from two test-particle simulations in 
isotropic turbulence with a nominal turbulence 
strength B = 1. Furthermore, the turbulent wave 
vector was scaled so that either k n • £ n = or the 
orientation of k is isotropic (see Sec. l2.4[) . However, 
such has only marginal influence on the resulting 
transport parameters compared to the estimated 
errors. 

In contrast, comparison of the classic results by 
Giacalone & Jokipii [6| , where no turbulence renor- 
malization had been done, with the Padian results 
shows systematical deviation, as clearly exhibited 
by Fig. [T] Therefore, the main deviation results 
from the fact that, in the Padian code, the turbu- 
lent magnetic field components have been normal- 
ized so that B — 1 on average. 

Moreover, analytical results that have been de- 



rived using second-order quasi-linear theory 24, 2 



agree better with test-particle simulations in a tur- 
bulent field with the correct turbulence strength. 
Such can be understood from the well-known fact 
that, as a rough estimate from classic quasi-linear 
theory one has A|| oc (B / Bq)~ 2 , thus underlin- 
ing the important influence of the magnetic field 
strength on transport parameters. Here, Bq = 
B^e. z denotes the mean magnetic field, which is usu- 
ally assumed to be homogeneous. 

To conclude, using the conventional approach of 
superposing plane waves, it is not possible to create 
a strictly isotropic turbulent magnetic field struc- 
ture that obeys all physical constraints, which are 

(i) equal mean of all magnetic field components; 

(ii) isotropy of the wave vectors; and (iii) vanish- 
ing divergence of the magnetic field. Such mag- 
netic fields are widely implemented in test-particle 
Monte-Carlo simulations, which are used to obtain 

(i) scattering mean free paths of charged particles; 

(ii) field line diffusion coefficients. 

While the turbulent magnetic field strength plays 
an important role for the results, such does not 
seem to be the case for a non-zero magnetic field 



divergence and/or the isotropy of the wave vectors. 
Future work should explore the possibility of a tur- 
bulence approach that is sufficiently simple but is 
fully compatible with all physical boundary condi- 
tions. 
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